Magnon-assisted Andreev reflection in a ferromagnet-superconductor junction 
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We study subgap transport in a ferromagnet-superconductor junction at low temperature due to 
Andreev reflection. The mismatch of spin polarized current in the ferromagnet and spinless current 
in the superconductor results in an additional contact resistance which can be reduced by magnon 
emission in the ferromagnet. Using the s-f model and focusing on half-metallic ferromagnets, we 
calculate the corresponding nonlinear contribution to the current which is asymmetric with respect 
to the sign of the bias voltage and is related to the density of states of magnons in the ferromagnet. 
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Spin polarized transport is a subject of intense re- 
search, motivated largely by the possibility of developing 
a form of electronics which utilizes the spin polarization 
of carrierso. Ferromagnetic (F) metals have more carri- 
ers of one spin polarization (known as majority carriers) 
present at the Fermi energy Ep than of the inverse po- 
larization (minority carriers). Of particular interest are 
so called half-metals where the splitting T between the 
majority and minority conduction bands is greater than 
Ep measured from the bottom of the majority band. In 
such a material only majority carriers are present at the 
Fermi energy and electric current in it is completely spin 
polarized. An F-F junction with antiparallel spin polar- 
izations has a larger contact resistance than a junction 
with parallel polarizations due to a mismatch of spin cur- 
rents at the interface, leading to tunneling magnetoresis- 
tance in F-F junctions™!] and giant magnetoresistance in 
multilayer structuresQtl Spin current mismatch may also 
affect the conductance of a ferromagnet-superconductor 
junctionErtl At low temperatures and small bias volt- 
age, current flows through the interface due to Andreev 
rcflcctionEj whereby particles in the ferromagnetic region 
with excitation energies e smaller than the superconduct- 
ing gap energy A are reflected from the interface as holes. 
Since subgap transport in the superconductor (S) is me- 
diated by spinless Cooper pairs the spin current is zero 
in the superconductor in contrast to the ferromagnet. 

Spin relaxation processesollij, such as spin-orbit scat- 
tering at impurities or magnon emission, can reduce the 
spin current mismatch. In a given junction, spin-orbit 
scattering (which is an elastic process) wojtJd reduce the 
value of the additional contact resistanceB'tHl whereas the 
inelastic process of magnon emission would manifest itself 
as a modification of the form of the I(V) characteristics. 
Nonlinear I(V) characteristics due to magnon-assisted 
tunneling between two feppmagnets have already been 
studied both theoretically^ and experimentallyEj with a 
view to relate the second derivative of the current to the 
density of states of magnons as d 2 I/dV 2 cx fl(eV). 

In the present work we analyze the role of magnon emis- 
sion in the formation of the subgap I(V) characteristics, 
eV < A, of a ferromagnet-superconductor junction with 
emphasis given to half- metallic ferromagnets. Microscop- 



ically, the relevant process involves the transfer of a sin- 
glet pair of electrons from the superconductor to the fer- 
romagnet where one of them forms an intermediate state 
at an energy below the bottom of the conduction band 
for minority electrons. Then, this latter electron can emit 
a magnon in the ferromagnet which will carry away ex- 
cess spin allowing the electron to incorporate_jtself into 
the majority conduction band. The s-f modet§E3 is em- 
ployed to describe a ferromagnet assuming that electrical 
conduction and magnetism are caused by different groups 
of electrons which interact via s-f exchange. 

To anticipate a little, due to magnon-assisted Andreev 
reflection at zero temperature, the I(V) characteristics 
of an F-S junction for eV < A acquire an additional 
contribution which is specified by the differential of the 
nonlinear differential conductance, 
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Here F(x) — x/(l — x ) , u> m is the magnon bandwidth 
and ujq is the anisotropy energy responsible for a gap 
in the magnon density of states JQ(w) is non-zero for 
uiq < ui < uj m and zero otherwise] t3. The factor of two 
in f2(2eV) arises because an elementary Andreev reflec- 
tion process involves a net transfer of charge 2e across 
the biased junction. The conductance Ga incorporates 
tunneling parameters relevant for the Andreev process 
and it is defined explicitly later. W is the conduction 
band width of electrons, v is the volume of a unit cell, 
and r = 2JS is the conduction electron band splitting 
where S is the spin of the local moments, and J is the 
s-f exchange constant. 

The nonlinear contribution to the current, Eq. ([!]), is 
asymmetric with respect to bias voltage V. It is zero for 
negative bias and finite for V > Uo/(2e). This feature 
can be expained using the sketch in Fig 1 which illus- 
trates the tunneling process between an S electrode on 
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the left hand side and an F electrode on the right for (a) 
V > and (b) V < 0. We have adopted the convention 
that positive (negative) voltage results in a Fermi energy 
Ep in the ferromagnet that is lower (higher) by energy 
\eV\ than the Fermi energy E$ in the superconductor. 
For V > 0, Fig 1(a), Andreev reflection results in the 
injection of both a majority (spin 'up') and a minority 
(spin 'down') electron into the ferromagnet. One of these 
electrons has an energy e above the Fermi energy in the 
superconductor E$, the other has an energy e below E$- 
With respect to the Fermi energy in the ferromagnet Ep 
these energies are eV + e and eV — e, respectively, which 
are both above Ep because of the need to move into un- 
occupied states in the ferromagnet, e < eV. At zero 
temperature T = the core spins are all aligned in the 
majority, up direction and, if Ep <C T, then only spin 
up conduction electrons are present at Ep in the ferro- 
magnet. A dynamic process which allows a spin down 
electron to enter the ferromagnet is shown schematically 
in Fig 1(a). The spin down electron tunnels from the su- 
perconductor into a virtual, intermediate spin down state 
above Ep , then it emits a magnon and incorporates itself 
into an empty state in the majority conduction band. In 
Fig 1(a) the magnon is depicted as a flip in the spin of 
one of the localized magnetic moments. 

On the other hand, for V < 0, Fig 1(b), Andreev re- 
flection would result in the injection of a spin up and a 
spin down electron from the ferromagnet into the super- 
conductor with energies e above and below Eg. With 
respect to the Fermi energy in the ferromagnet Ep these 
energies are — \eV\ +e and — \eV\ — e, respectively, which 
are both below Ep because of the need to have initially 
occupied states in the ferromagnet, e < \eV\. As stated 
above, there are no spin down states at these energies 
near Ep. Since a spin up electron cannot emit a magnon 
(left side of Fig 1(b)), due to conservation of total spin 
in the exchange interaction, the only possibility would be 
that a spin up electron in the ferromagnet would absorb 
a magnon to form an intermediate spin down state before 
tunneling into the superconductor. However there are no 
thermally excited magnons at T = in the initial state of 
the ferromagnet (right side of Fig 1(b)) so it is impossible 
for magnon-assisted Andreev reflection to contribute to 
current formation in negatively biased junctions. 

To calculate the current across an F-S junction, we use 
the transfer HamiltaMAn method and the nonequilibrium 
Keldysh technique£§ll!a to express the current in terms of 
the tunneling density of states in the metallic ferromag- 
net which will be determined using the s-f model. We will 
generalize the approach of Ref. |l9| developed for normal- 
superconductor junctions to the case of an F-S junction. 
We consider a point contact geometry which consists of 
a narrow constriction between S and F electrodes and we 
use the total Hamiltonian of the system in the form of 



H — Hg + Hp + Hp, 
H T = y~] \tc\ a c F , a + 1* c F a c s ., 



in all calculations paying attention to the spin index 
a — {T: I}- 111 the tunneling Hamiltonian Hp, t describes 
the transfer of an electron from S to F. We assume that t 
is independent of energy and that spin is conserved upon 
electron transfer across the interface. 

Assuming that the voltage drop V occurs across 
the point contact and both electrodes are almost in 
equilibrium, for a given distribution and density of 
states in the reservoirs, the current can be calcu- 
lated as the mean value of the current operator, / = 

(ie/h) X) a (* ( C s,a c -F>) ~ ** ( c pa c S.a))- The averages in 
this equation may be expressed in terms of nonequilib- 
rium Keldysh Greens functions G" 1 and it is convenient 
for the case of a superconducting electrode to use the 
2x2 Nambu representationEl After some algebra we 
arrive at the following expression for the current 
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This expression describes the current due to Andreev 
reflection of electrons at the F-S interface. The gap- 
full spectrum of quasi-particles in the superconductor is 
taken into account in the off-diagonal element of the re- 
tarded Greens function Gg 12 and the energy e is mea- 
sured with respect to the Fermi energy Es in the super- 
conductor. 

In Eq. (||), p a {E) is the tunneling density of states 
of spin state a in the ferromagnet. For T ^> {Ep ,eV} 
and in the absence of spin waves, Pi(E) is zero in a half 
metal at energies close to the chemical potentials of the 
electrodes and I(V) = 0. However, as explained above, 
Pl(E) is nonzero for E > Ep + ui because of the dy- 
namical process of magnon-assisted tunneling. The tun- 
nelina-deiisity of states pi (E) is determined using the s-f 
mode]li§E3 with the following Hamiltonian 



Hp = H o + Hj + H ex , 
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where Hq deals with conduction band electrons, Hj is 
the Heisenberg Hamiltonian of the localized moments and 
H ex takes into account the s-f intra-atomic exchange. We 
will consider weak, ferromagnetic s-f coupling < J S < 
W, thus excluding the case corresponding to double ex- 
change JS > W. 

The conduction electron single particle Greens func- 
tion and energy spectrum have been studied in the 
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limit p£ zero temperature-both for zero conduction band 
fillingjlj and finite fillingEj. In the regime {eV,LOo} <C 
Ep <C r <C W the minority electron density of states is 
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where the energy E is measured from the bottom of the 
spin up conduction band and L d is the volume of the 
system. A similar expression, with 2S — 1, was also 
obtained for itinerant ferromagnets using the Hubbard 
modelEa. The energy in the denominator of Eq. (||), 
f-k'+q — E + T, is related to the inverse lifetime of an 
electron in the virtual state, the delta function in the 
numerator accounts for energy conservation of the en- 
tire process and the factor [1 — f(k')] indicates that this 
process is a dissipative one and only involves unoccupied 
final (spin up) states. Therefore the minority density 
of states is finite for energies E > Ep + loq only. At 
T = 0, in our regime of interest, spin waves do not af- 
fect the majority density of states P]{E) ps N(E), where 
N(E) is the bare electron density of states arising from 
the Hamiltonian Hq. We have neglected interactions be- 
tween conduction electrons so that the dominant effect of 
finite filling is the requirement to obey Fermi statistics. 

The dispersion of magnons in a ferromagnet is 
quadratic at low frequencies, ui q = Dq 2 + loq where 
D = 2JS 2 a 2 and a is the lattice constant. For sim- 
plicity, we shall also use the Debye approximation for a 
broader range of energies of the magnon spectrum by in- 
troducing a cutoff of the parabolic magnon dispersion at 
uj m ps D(6tt 2 /v) 2 / 3 since, for A <C uj m , such a simplifi- 
cation does not influence the final results. We also use 
a parabolic approximation for free electrons in a conduc- 
tion band and find that, for a three dimensional ferro- 
magnetic metal and E ss Ep , 
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We use Eq. (Q) to calculate the current due to Andreev 
reflection for low bias voltage eV < A, taking into ac- 
count the lowest order in the tunneling matrix elements, 
\t\ 4 . In a BCS superconductor the off-diagonal Greens 
function is Gf. 12 (e) » ttN{E s )/^1 - (e/A) 2 for e < A, 
where N{Eg) is the density of states in the normal state 
near the Fermi energy. Note that if the ferromagnet were 
to be replaced by a normal metal where the density of 
states at Ep is N{E F ) for both spin channels, the cur- 
rent at low bias would be / = VG a with conductance 
G A m (e 2 /h)16n 4 \tf N 2 {E F )N 2 (Es). As we assume 
eV < E F , p T fa N(E F ), and Eq. (§) for an F-S junc- 
tion can be simplified, by shifting the energy variable, 
into 
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we find for eV < A <C 
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The I(V) characteristics obtained yield the relation be- 
tween the magnon density of states and d 2 I/dV 2 shown 
in Eq. (|j). For small voltages eV <C A, the current may 
be simplified further as 
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and d 2 I/dV 2 cx fi(2eV). The factor 6(2eV-oj ) empha- 
sises the absence of a contribution to the current from 
magnon- assisted Andreev reflection for negative bias. At 
zero temperature, the /(V) characteristics are asymmet- 
ric, which can be used to extract magnon-assisted An- 
dreev reflection in an F-S junction involving a ferromag- 
net with electrons of both polarizations present at Ep. 
Note that at T > u>o, the absorption of thermal magnons 
would allow for a similar contribution to the current at 
V < 0. A graphic representation of the result of Eq. (|l|) 
is given in Fig 2, in comparison with a similar analysis for 
an F-F junction in an antiparallel configuration (inset). 
In particular, we illustrate how the relation between the 
conduction electron bandwidth and the exchange split- 
ting energy, W/T, may affect the result, since for W ^> F 
the contribution of large energy magnons to the lifting of 
spin current mismatch is not so efficient as for W <C F. 

To summarise, spin current mismatch in an F-S junc- 
tion at subgap voltages may be lifted by magnon-assisted 
tunneling resulting in nonlinear /(V) characteristics re- 
lated to the magnon density of states. Such features are 
typical not only for junctions between ferromagnets and 
BCS superconductors, but also for HiTc materials, since 
subgap transport through any spin singlet superconduc- 
tor would be affected by spin current mismatch at the 
F-S interface. For a matching pair of materials, with the 
Curie temperature of the ferromagnet similar to the tran- 
sition temperature of the HiTc material, the Andreev re- 
flection process could be used to probe the whole magnon 
spectrum. 
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FIG. 1. Schematic of the tunneling process between 
a superconducting electrode on the left hand side and a 
ferromagnetic electrode on the right for (a) V > and 
(b) V < 0. For (a) V > a down spin electron tunneling 
into the ferromagnet may emit a magnon and incorporate 
itself into the majority conduction band. For (b) V < 
no spin flip process is possible at T = because in the 
initial state (right) there are no thermal magnons for an 
up spin electron to absorb. See text for details. 
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FIG. 2. Typical form of d 2 I/dV 2 for an F-S junction 
with various choices of material parameters. From top 
to bottom, Y/W — 5.0 (dot-dashed), 1.0 (dashed), and 
0.2 (solid). As an example, we choose luq = 0.2A and 
u) m = 4. OA. For comparison, inset shows d 2 I/dV 2 for an 
F-F junction with antiparallcl spin polarizations. 



